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Abstract 

An energy eigenfunction in a classically chaotic system is known to have 
spatial correlations which (in the limit of small K) are governed by a micro- 
canonical distribution in the classical phase space. This result is valid, how- 
ever, only over coordinate distances which are small compared to any relevant 
classical distance scales (such as the cyclotron radius for a charged particle in 
a magnetic field). We derive a modified formula for the correlation function 
in the regime of large separation. This then permits a complete description, 
over all length scales, of the statistical properties of chaotic eigenfunctions in 
the h ^ limit. Applications to quantum dots are briefly discussed. 
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In a hamiltonian system which exhibits classical chaos throughout the accessible phase 
space, the quantum energy eigenvalues and eigenfunctions are known to have certain uni- 
versal properties in the limit of small H (which, in practice, is achieved at sufficiently high 
energy) [|]||. We will be interested in the energy eigenfunctions; these can be characterized 
as random variables with a gaussian probability distribution of the form p|-p 



P(ip\E) oc exp 



| J dki J dfq 2 V>*(q 2 )#(q2,qiW(qi) 



(1) 



where / is the number of degrees of freedom. If the system is time-reversal invariant, 
the eigenfunctions are real and (3 = 1; if it is not, they are complex and (3 = 2. The 
measure corresponding to eq. ([I]) is the standard one of euclidean quantum field theory: 
Dip = n q #(q) for /5 = 1 and Vip = n q dRei/j(q) dim ^(q) for (3 = 2. P(^\E)V^ 
represents the probability that the actual eigenfunction ip a (<i) for E = E a (a particular 
energy eigenvalue) is between ^>(q) and ip(q) + dip(q) for all coordinates q. The kernel 
^(q2> qil-E 1 ) is the functional inverse of the two-point correlation function 

C(q 2 , qi |£) = J MthWMPWEW ■ (2) 

The explicit formula for C(q 2 , qjJ-E 1 ) which was originally suggested by Berry |3|] assumes a 
microcanonical probability density in the classical phase space, 

C(q 2 ,qi|£) = ^E)j jP<w»H(E - H w (p,q)) , (3) 

where q = |(qi + q 2 ), Hw(p, q) is the classical hamiltonian (more specifically, it is the Weyl 
symbol of the hamiltonian operator), and p(E) is the semiclassical density of states, 

However, eq. (|3|) only applies when the separation |q 2 — qi| is sufficiently small M. Consider, 
for example, the case H = p 2 /2m + V^(q); we would not expect the correlations in an 
eigenfunction at two points qi and q 2 to depend only on the value of the potential at q if 
V(q) has a significantly different value than either V^(qi) or V^(q 2 ). Our goal, then, is to 
find the correct formula for C(q 2 , qi\E) when |q 2 — qi| is large. 



One way to motivate eq. (Rl) is to consider the energy Green's function flP2"|,|l"3 



r(n „ _>^ ^(q 2 )C(qi) f ,s 

G(q 2 , qi |E)=E E _ (5) 



where e — > + . We then have 



E ^(q 2 )C(qi W£ ~E a ) = ^- [G(q 1; q 2 \EY - G(q 2 , q^E)] , (6) 

2m 



and the exact density of states is 
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p(E) = £ 5(E - E a ) (7) 

a 

While the exact Green's function clearly has singularities whenever E = E a , its leading 
approximation G(q 2 , qi|-E) in the small-/;, limit is a smooth function of its arguments. Given 
this, eqs. (^,^,0) make it natural to expect that 

C(q 2 , qi |£) = 2 J p^^E)* -G^^E)] . (9) 

This formula can be derived explicitly || in the theory of disordered metals, where a white- 
noise random potential is added to the hamiltonian. In this case, eq. @ holds with G 
standing for the Green's function averaged over the random potential. In the limit that 
the strength of the potential is large (which corresponds to the % — > limit for chaotic 
systems), eq. (P for the eigenfunction probability can also be derived explicitly |5H8l,|10 



Corrections for finite potential strength (that is, finite Ti) can also be computed W% , and 
are complimentary to the calculations done here. 



Let us now briefly recall the construction of G(q 2 , qil-E 1 ) in the Ti — > limit [jl2|,[L3 . The 
energy Green's function is related to the propagator (q 2 |e~* H */^|qi) via 

1 r°° 

G{q 2 , qi \E) = - dt e^+^qsle-^lqx) . (10) 
in Jo 

The propagator can be written as 

(q 2 |e-^| qi ) = / j^j e^-^^e-^V(P,q) , (H) 

where again q = |( qi + q 2 ), and Avi/(p,q) denotes the Weyl symbol of the operator A; 
in fact, eq. (0) is simply the Fourier transform of the definition of the Weyl symbol. For 
sufficiently small times, 

{e- iHt/h ) w {p, q) - e - iHw{ ™ )t/h . (12) 



Inserting eqs. (|ll~l , |l2"l) into eq. fllOl), and performing the time integral, we get the desired 
approximation, 

G(q 2 , qi |£) = / -^-j e *P-fo»-iO/ R —1 . (13) 

^H2,hii ) j E-H w {p,q)+ie K ' 

The derivation of this formula is flawed, however, since we integrated over all positive times 
even though eq. (|1^) is valid only for short times. A more careful analysis shows that eq. (|13|) 
is valid in the limit of small h, provided that |q 2 — qi| is small enough so that the shortest 
classical path q(t) connecting qi to q 2 with energy E is well approximated by a linear 
function of time. Note that this criterion is purely classical; |q 2 — qi| can be sufficiently 
small even if it is large compared with the quantum wavelength h/p, where p = |p| is the 
magnitude of the classical momentum when q = q. With this caveat, eq. (|T3|), when inserted 
into eq. (P), immediately yields eq. (|3]). 
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As this derivation shows, however, eq. @ is not valid if \q% — qi| is too large. In this 
case, we must use a different semiclassical formula for G(q2, q.i\E), making a stationary 
phase approximation both in the Feynman path integral representation of the propagator, 
and in the time integral of eq. (|T0|) . The well-known result is [12[ 



G(q 2 ,qi|£) = .} Mf _ 1)/2 E iD^e^-^ 2 . (14) 



ih[2~;i,yf 1 



Here the sum is over all classical paths connecting q x to q2 with energy E and action 



/ p-dq. (15) 



The index z/ p counts the number of classical focal points along the path, and the determinant 
D p of second derivatives of S p is given by 



D p = det 



( a 2 s p a 2 s p 

9q2<9qi dEdqi 

d 2 S p d 2 S p 

\ dc^dE dE 2 



(16) 



It is eq. (|14D which should be used in eq. (Q) when the shortest classical path from qx to q2 
is not (approximately) a linear function of time. 

If the system is time- reversal invariant, then the paths from q2 to qi are have the same 
set of values of S p , u p , and D p as the paths from qi to q2. In this case, we have 

C{q>ME) = p (E)(27r ft)(/ + i)/2 £ l^ P r /2 cos[S P A-(2^ P + /-lW4] (17) 

instead of eq. ([|) when | q 2 — qi | is large. 

To illustrate the differences between eq. ([[7]) and eq. @), let us examine a few special 
cases. First, we consider an /-dimensional billiard in which the straight-line path from qi 
to q2 is not blocked. In the interior of the billiard, H = p 2 /2m, and eq. @ yields |J 

C^ME) = V-^UI2) J ( l'-J^ l2 , (18) 

where V is the /-dimensional volume of the billiard, hk = (2mE) 1 ^ 2 , L = |q 2 — qi| is the 
length of the straight-line path, T(x) is the Euler gamma function, and J v (x) is an ordinary 
Bessel function. In the case at hand, eq. (O) is valid even for large L. The reason is that 
the straight-line path is a linear function of time, and this is the condition needed for the 
validity of eq. (gj. Turning to eq. (0), we note that p(E) = yV/{ATi)f' 2 Y{f/2)E, and 
that the straight-line path has S p /h = kL and \D V \ = m 2 {2mE)^~^/ 2 / 'i/' 1 . If q x and 
q2 are both far from any of the billiard's walls, the straight-line path makes the dominant 
contribution, and we find 

C[&ME)-V r U/2) 7r i/2 (fcL/2 )(/-i)/2 ( 19 ) 

with z/p = 0. Eq. (|19|) is equivalent to eq. ( |18D when kL is large, since in this regime 
the asymptotic form of the Bessel function can be invoked. Thus, in the present case, both 
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eq. (Q) and eq. (|T7|) are valid for large L, and eq. (|) is valid for small L (less than a quantum 
wavelength, 2ir/k) as well. 

On the other hand, if we consider a billiard in which the straight-line path from q x to q2 
is blocked by an obstacle, such as in fig. (1), eq. (|18D is not correct, and we must use eq. (|17|). 
The shortest classical path connecting qi to q2 makes the dominant contribution; this is 
given by eq. (|19|), provided we take L to be the length of the path, and set v v equal to twice 
the number of bounces. In most realistic cases of this type, there will be many other paths 
with more bounces that are not very much longer; these will all contribute to C(q2, qi|D) as 
well. However, if we attempt to verify eq. (|TjJ) numerically for a particular eigenfunction of 
a particular billiard, it is necessary to average over a range of qi and q2 in order to reduce 
the variance in C(q2,qi|£') which is expected from the probability distribution (|l|) PfllTTj . 
If this averaging is carried out with the length L of the shortest classical path held fixed, 
the other contributing paths will in general have lengths that vary. If this variation is large 
on the scale of the quantum wavelength 27r/fc, then the net contribution of all these other 
paths to the averaged C(q2,qi|-D) should be small, rendering eq. (|T9|) a valid formula for 
the averaged correlation function. 

As another example we consider a billiard with the addition of an isotropic harmonic 
potential V(q) = |mco> 2 q 2 . This case is of some physical interest for / = 2; wave-function 
correlations in quantum dots have been studied assuming that the dot is well-modeled by 
a two-dimensional billiard, but in fact there is also a smooth confining potential, often 
approximated as harmonic. For simplicity, we give only the leading corrections in the limit 
of a weak potential. The shortest classical path (assuming it is not blocked) has action 



S p /h = kL 



L 2 - 3d 2 
1 + 12R 2 + 



(20) 



where hk = {2mE) l l 2 and L = | q 2 — Qi I as before, and we have introduced d 2 = q 2 + q 2 , 
and R = (2E /mw 2 ) 1/2 i R is the maximum distance from the origin which can be reached 
with energy E. Also, we find 



I -Dp I — | -Dp 1^=0 



{f - 1)L 2 - (/ - 3)d 2 
AR 2 



(21) 



In all but two dimensions, the correction to \D P \ is dominated by the correction to p(E), 
which is 0(V 2 /f / R 2 ), where V 1 ^ is the linear size of the billiard. However, when / = 2, 
p{E) is independent of E, and it is not changed by the presence of a weak potential. In this 
case, eq. fl2T|) represents the dominant correction to the amplitude of C(q2,qi|D). 

As a final example with importance for quantum dots, we consider a particle with charge 
e in a two-dimensional billiard with a uniform, perpendicular magnetic field B. In the 
billiard interior we have H = (p — eA) 2 /2m, and we will work in the gauge in which the 
vector potential is A = |B x q. This system is not time-reversal invariant, and so we must 
use eqs. (|9]) and flU]) rather than eq. (|iT|). We again consider points qi and q 2 which are 
far from the billiard's walls. The shortest classical path is then a circular arc with length £, 
related to the separation L = |q 2 — qi| and classical cyclotron radius R = {2mE) l l 2 /\eB\ 
via £ = 2Rsin~ 1 (L/2R). The action for this path can be divided into a geometric part and 
a gauge- dependent part, S p = S gcom + S gange . The geometric part is 
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where again hk = (2mE) 1 ^ 2 , and A = \Rt — \R 2 sm(£/R) is the area enclosed by the 
circular arc and the straight line connecting qx to q 2 . The gauge-dependent part is energy 
independent, and changes sign when q! and q 2 are exchanged. For our gauge choice, 

•S'gauge = |eB-(qi x q 2 ) . (23) 

If we make a gauge transformation 

A(q)^A(q) + V$(q), (24) 

where $(q) is any smooth function, then 

Sgauge Sgaugc + e$(q 2 ) - e$(qi) . (25) 

The determinant \D P \, on the other hand, is gauge invariant, 

^^Mi^-Tm) ■ (26) 

Again there is no correction to p(E) in two dimensions. Keeping only the contribution of 
this path, we find from eqs. (||) and QT4]) that 

C (q2 , qi | B ) = V 1 ^W ,,^^ • (27) 

where V is the area of the billiard. Under the gauge transformation (p4|), eq. (|25|) implies 

C(q 2 , qi \E) - e + ^^-^/ n C(q 2 , q,\E) . (28) 
That this is correct can be seen by recalling that a wave function ^(q) transforms as 

V>(q) -> e +ie * (q)/ V(q) (29) 



under (E3), and that C(q 2 , qi|-E a ) is the expected value of ^«(q2)V'a( c li)- On the other 
hand, eq. (§) implies that 

C(q 2 , qi |£;) - e +^-^ v ^C(q 2 , qi |£) , (30) 

which again illustrates the fact that eq. (|3]) is valid only when |q 2 — qi| is sufficiently small. 

Finally, we note that our expression for C(q 2 ,qi|-E) is needed to resolve a discrepancy 
between two different formulas in the literature for the expected values of off-diagonal matrix 
elements (in the energy-eigenstate basis) of simple, ^.-independent operators in classically 
chaotic systems. This will be the subject of a separate paper. 



ACKNOWLEDGMENTS 

This work was supported in part by NSF Grant PHY-97-22022. 



6 



REFERENCES 



[1] Les Houches LII, Chaos and Quantum Physics, ed. M.-J. Giannoni, A. Voros, and J. 

Zinn- Justin (North-Holland, Amsterdam, 1991). 
[2] Quantum Chaos, ed. G. Casati and B. Chirikov (Cambridge Univ. Press, New York, 

1995). 

[3] M. V. Berry, J. Phys. A 10, 2083 (1977). 

[4] P. O'Connor, J. Gehlen, and E. J. Heller, Phys. Rev. Lett. 58, 1296 (1987). 
[5] V. N. Prigodin, K. B. Efetov, and S. Idia Phys. Rev. Lett. 71, 1230 (1993). 
[6] V. N. Prigodin, B. L. Alshuler, K. B. Efetov, and S. Idia Phys. Rev. Lett. 72, 546 
(1994). 

[7] V. N. Prigodin, Phys. Rev. Lett. 74, 1566 (1995). 

[8] V. N. Prigodin, N. Taniguchi, A. Kudrolli, V. Kidambi, and S. Sridhar, Phys. Rev. Lett. 
75, 2392 (1995). (|cond-mat/9504089|) 

[9] Y. Alhassid and C. H. Lewenkopf, Phys. Rev. Lett. 75, 3922 (1995). ( |cond-mat/9510023| ) 
[10] M. Srednicki, Phys. Rev. E 54, 954 (1996). (|cond-mat/9512115|) 
[11] M. Srednicki and F. Stiernelof, J. Phys. A 29, 5817 (1996). ( |chao-dyn/96030l2| ) 
[12] M. C. Gutzwiller, J. Math. Phys. 8, 1979 (1867). 
[13] M. V. Berry and K. E. Mount, Rep. Prog. Phys. 35, 315 (1972). 
[14] Ya. M. Blanter and A. D. Mirlin, Phys. Rev. E 55, 6514 (1997). (|cond-mat/9604139|) 



7 



FIGURES 




FIG. 1. A portion of a Sinai billiard in which the circular scatterer is off center; the path 
shown is the shortest classical path between points one and two. 
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